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This paper reviews the Dalitz plot analysis for the decays B'\B'^) pTv ^ Tr+Tr^Tr". We discuss 
what can be learned about the ten parameters in this analysis from untagged and from tagged time- 
integrated data. We find that, with the important exception of the interesting CP violating quantity 
a, the parameters can be determined from this data sample - and, hence, they can be measured 
at CLEO as well as at the asymmetric B factories. This suggests that the extraction of a from 
the time-dependent data sample can be accomplished with a smaller data sample (and, therefore, 
sooner) than would be required if all ten parameters were to be obtained from that time-dependent 
data sample alone. We also explore bounds on the shift of the true angle a from the angle measured 

, from charged p final states alone. These may be obtained prior to measurements of the parameters 

• describing the neutral p channel, which are expected to be small. 
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I. INTRODUCTION 



The B factories at SLAC and KEK are now collecting data and expect to produce measurements of the CP 
asymmetry in the mode B — > J /ipKs within one year, measuring sin 2/3. Preliminary results on this mode from CDF 
already indicate that it is unlikely that a discrepancy with the Standard Model will be found from this result 
alone. This means that tests of the Standard Model mechanism for CP violation will rely upon our ability to measure 
further CP- violating parameters, such as the angle a = 7r — /3 — 7of the Unitarity triangle, with sufficient accuracy to 
^-H ' be sensitive to Standard Model relationships. This will require that we master the removal of theoretical uncertainties 
, due to penguin diagrams in at least one of the available channels. So far there are two sets of candidate decay modes, 
' TTTT and pTT 1^], for which analyses to extract a using isospin relationships have been suggested. The first suffers 
\ from relatively small branching ratios [Q and from the experimental difhculty of measuring the 7r°7r° branching ratio. 
(~| • Some of the modes for the second case have recently been observed at CLEO § . These results are encouraging. They 
rS i' are close to model-dependent predictions and, hence, appear to reinforce estimates that the analysis suggested by 
JL Snyder and Quinn will require of order 180 fb~^, or six years of running at BABAR design luminosity, to complete 

, This paper revisits that analysis and reviews what intermediate steps can be made with earlier data samples. In 
IL" ■ particular we stress that many of the parameters relevant to the extraction of a can be predetermined from untagged, 
and from tagged but time-integrated data samples. The benefit of this is that, once this is done, only the CP- violating 
parameter a remains to be fit to the full tagged and time-dependent Dalitz plot. Presumably it will then be possible 
^ I to perform the pir analysis with smaller data samples than those needed if the full set of ten parameters is fit to the 
tagged, time-dependent sample alone (as has been done to date in Monte Carlo studies of this mode). We also discuss 
bounds on the shift in a from penguin contributions. These can be obtained, by methods similar to those previously 
suggested for the tttt modes , even if the p^Tr" amplitude is too small to be measured directly. We will also discuss 
a bound which applies only in the pir case, having no parallel in the tttt case. 

This paper presents a purely theoretical discussion, with no simulations and no attempt to address issues of back- 
grounds or of other modes that may contribute to the 7r+7r~7r° Dalitz plot ||l^. Certainly these issues will be 
important. The larger untagged data sample will also be the first place to explore the issues of further contributions 
to the Dalitz plot. The approach of fixing as many parameters as possible from the untagged and from the tagged 
but time-integrated data samples, combined with improvements in machine luminosity, such as those already under 
study for PEPII, may make the extraction of a reliable value for a a reality on a somewhat faster time scale than 
suggested by the estimates based on fits to tagged data only. 



*Permanent address: Instituto Superior de Engenharia de Lisboa, Rua Conselheiro Emfdio Navarro, 1900 Lisboa, Portugal. 
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II. NOTATION 



A. Isospin decomposition 

We are interested in the decays from and 5° into pn final states. The decay amphtudes can be classified 
according to isospin: {B~^, B^} form an isospin doublet; the final state pn can have isospin // — 0, If = 1, and If = 2 
components. In general, the final state with If ~ can only be reached with operators having A/ = 1/2; the final 
state with If = 1 can be reached with operators having A/ = 1/2 or A/ = 3/2; and the final state with If = 2 can 
be reached with operators having A/ — 3/2 or A/ — 5/2. We denote the isospin amphtudes by A^ijf Thus, 



a+o ^a{B+ ^ p+t:°) = ^^^^3/2,2 - ^-^^3/2,1 + ^^1/2,1 



ao+ = a{B+ ^ p\+) 



a+- ^a(B° ^ p+TT' 



a (B° ^ p^7r^ 



1 /3 



A, 



1 1 



2 V 2-^3/2.2 +2^^/2,1 



1 



^^3/2,2 + 2^3/2,1 + 2^1/2,1 + ;^^V2,0 



1 



= ^3 



2^^i3/2,2 - 2^3/2,1 - 2^1/2.1 + 71^1/2,0 



1 



^.^3/2,2 - ^A/2.0 



-^AV2,2 



1 



2^/3 
1 



5/2,2 



(1) 



Similar relations hold for the decay amplitudes, which we denote by a_o, oo-j o hi '^H > ^^'^ '^00 1 respectively. 

Note that our notation here is that an is the amplitude for the B^ to decay to a p of charge i and a tt of charge j. 



The CP relationships are thus CP{aij) — 



The isospin components are Aaij^ ■ 



In the Standard Model, there are tree- level amplitudes, gluonic penguin amplitudes, electroweak penguin ampli- 
tudes, and final state rescattering effects. The tree level b uud decays have both A/ = 1/2 and A/ = 3/2 
components. In contrast, the gluonic b d penguins are pure AI = 1/2, because the gluon is pure / = 0. Therefore, 
the isospin amplitudes ^13/2,1 and ^3/2,2 only receive contributions (and, thus, weak phases) from the tree-level dia- 
grams. There are no diagrammatic AI = 5/2 contributions at this order; such effects arise only from electromagnetic 
corrections to the weak-decay diagrams. 

A priori, there is no hierarchy among the AI =1/2 and AI = 3/2 isospin amplitudes. However, the combination of 
isospin amplitudes involved in the decay B'^ — > p^n'^ is generally argued to be suppressed because the tree-level and 
gluonic penguin diagrams can only contribute to B^ —> p^w'^ through a color-suppressed recombination of the quarks 
in the final state. This argument is not theoretically rigorous because B'^ p^ir'^ may be fed from other topologies 
through strong final state rescattering. Eventually experiment will tell us whether these effects are important. 

In what follows we will neglect two contributions which are suppressed in the SM. The first contribution arises from 
the electroweak penguin diagrams, which have the same weak phase structure as the QCD penguin diagrams, but 
contribute to both AI —1/2 and AI = 3/2 amplitudes. As a result, the electroweak penguin contributions are not 
removed by the isospin-based analyses. However, they are expected to be very small in these channels [pp]]. The 
second contribution is due to a possible AI — 5/2 isospin component, included within square brackets in Eqs. (|l|). In 
the SM, this component comes from electromagnetic rescattering effects and, thus, it is suppressed by a ^ 1/127.0 
Both effects could become relevant for 5° P^tt", should this amplitude turn out to be very small. 

Henceforth, we will only consider the tree-level and gluonic penguin diagrams. These give contributions with weak 
phases 



\y:,yt<i\ 



— g»7 — _g-i(/3+a) 



(2) 



'^Notice, however, that this effect is important in K - 
amplitudes, ReA2/Re^o ~ 1/22, encoded into the A/ 



TTTT decays because there one has a strong hierarchy among the decay 
1/2 rule llMl. 
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The first of these expressions is the weak phase of the tree-level contributions. Here we follow H and use the unitarity 
relationship to rewrite the penguin contribution as a dominant term proportional to the second of the CKM factors 
in Eq. plus a sub-dominant term proportional to the first CKM factor in Eq. In what follows we always 

subsume this second term within the amplitudes we refer to as AI —1/2 tree amplitudes. Since we do not calculate 
these quantities, but rather discuss extracting them from fits to experiment, this makes no difference to our analysis. 
Note, however, that if the amplitudes extracted in this way are to be compared to those calculated in any given model, 
then the penguin contributions to our so-called "tree" amplitudes must be taken into account. 

The interference between the two amplitude contributions in Eq. (H) depends on the weak phase a. However, we 
will show that, as is usual with direct CP-violation effects, any sensitivity to a is masked by an unknown coefficient 
with large theoretical uncertainty. Another weak phase arises from the interference between the i?° — mixing, 
q/p — exp[— 2i(/3 — 9d)], and the tree-level diagrams. For example, 

P ^3/2,2 

where dd parametrizes a possible new physics contribution to the phase in 5° — mixing. In the SM, 9d = and 
the two phases coincide; other models have 9d ^ and a difference arises. The phase probed by the Snyder-Quinn 
method is a + 6d . 

B. B° — > TT^Tr^Tr" decay amplitudes 
The /9-mediated B^ — > Tr+Tr^Tr" decay amplitudes may be written as 

where f± and /o are Breit-Wigner functions representing and p", respectively, and also include the cos 9 angular 
dependences of the hclicity zero p decays. The crucial observation made by Snyder and Quinn is that the Breit-Wigner 
functions contain CP-even phases, and that the interference between the different Breit-Wigner shapes across the 
Dalitz plot provides experimental sensitivity to the weak and strong phases contained in Eqs. (Q). There is some 
systematic uncertainty associated with the specific choice made for the shape of the functions f± and /q. This 
uncertainty affects primarily the corners of the Dalitz plot where the tails of two such functions interfere; these are 
the regions of the Dalitz plot from which one extracts the interference between two distinct decay amplitudes. 
It will prove convenient to rewrite Eqs. (||) as 

A = fcac + fdCld + /nfln, 

A^ fcac + fdad + fnan, (5) 



where 



and 



ac = a-i—+a-+, ad^a+--a-+, a„ = 2aoo, 

ac = a+-+a-+, ad = a+_-a_+, a„ = 2aoo- (7) 



Notice that the CP-conjugate of is not ad, but rather —ad- 

As discussed in the appendix, we parametrize the amplitudes Oc, ad, and a„ as 



^The term is sub-dominant in that it is a difference of up and charm quark contributions and, hence, it vanishes in the limit 
that these two quark masses are taken to be equal. 
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a, = Te-"" (1 - z - ro) , 
ad^Te-"" (zi+n), 

a„ = Te-^" (z + ro) , (8) 

where z and zi are CP-even, while tq and ri are CP-odd. The quantity Tzi contains the CP-even contributions to 
the final state with // = 1, summing the tree amplitude and the penguin amplitude multiplied by cosa. Tri {Ttq) 
contains the CP-odd part, given by the penguin contributions to the final state with If = 1 [If = 0), multiplied by 
I sin a. Similarly, the amplitudes contained in qAf/p can be written as 

^d^ = e^^'^'^Te'"' (1 - z -I- ro) , 
P 

^-ad = e'''-Te'^ (-^i+ri), 
P 

^a„ = e^^'-Te^^ (z - ro) . (9) 
P 

We have chosen to define strong phases so that T is a real positive quantity. Notice that Cc + a„ ~ Te"*", and 
q/p{ac + an) = e^^^-^Te*". Therefore, the imaginary (real) part of the ratio of these quantities, measures the sine 
(cosine) of the phase in Eq. (^) . 

If floo and aoo are indeed color suppressed, then z and ro will be smaller than 1, zi and ri. In that case, the 
CP- violating difference 

l„ |2 1^ |2 

= -2Rero + 2Re(zr5) (10) 



2T2 

will also be small. 

In connection with Eqs. (^ and (^, there is no advantage of one parametrization with respect to another. We 
could equally well have chosen to write the amplitudes in a new set of basis functions {/i, /2, /a}, given by 

f f+ + f- + 2/0 r f+ - f- r f+ + f- - fo -, X 

h — Q , J2 — ^ , J3 — ^ ■ (.-L-LJ 

In this basis, the amplitudes would become 

A = fi {ac + a„) + f2ad + h (flc - a„/2) 

fi (flc + an) + f2ad + h (a-c ~ a„/2) . (12) 

Eqs. (^ highlight the intermediate pn states; Eqs. (^) highlight the CP structure of the intermediate charged pn states, 
but still treat the intermediate p^n^ separately; and Eqs. ( |l2|) highlight the extraction of Oc -I- a„ = Te"'". Although 
each basis has its pedagogical advantages, they have no experimental significance. What one does experimentally is 
a maximum likelihood fit of the variables T, zi, ri, z, and ro to all the data in the Dalitz plots. 

C. Observables in the decay rates 

In the Bd system we have \q/p\ ^ 1 and AF <C F. Therefore, 

F[P" ^ /] + F[P^ ^ /] « \Af\' + \Af\^ (13) 
F[P° ^ /] - F[P0 ^ /] oc {\Aff - \Aff) cosAmi- 2Im (^^AfA)^ sin Ami. (14) 

One may rewrite these expressions using A/ = / (pA/ ) The three terms in Eqs. ( |T^ ) and ( p^ become 1 -I- |A/p, 
1 — |A/p, and ImA/, respectively. The untagged decay rate in Eq. ( p^ probes only the combination 1 -I- |Ayp, 
regardless of whether these measurements are time-dependent or time-integrated. Due to the anti-symmetric nature 



^Notice, however, that some authors use the opposite sign convention for q/p and, thus, for A/. 
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of the — pairs produced at the T(4s), tagged, time-integrated measurements performed at facilities working on 
this resonance can probe 1 + |A/p and 1 — |A/p, but not ImA^. 

If / is a CP eigenstate and the decay amphtudes Af and Af are defined in such a way that the rates and 
I Ay P include all the phase space integrations, then 1 + | Ay ^ is CP-even, while 1 — |A/p and ImA/ are CP-odd. Under 
these conditions, the ratio of Eq. ( |l^ to Eq. (|l3|) is the famous CP- violating asymmetry. It is sometimes stated that 
one can only probe ImA/ /(I + |A/|^. While this is true if one looks only at the CP- violating asymmetry, one can see 
from Eqs. ( [l^ ) and (|lj) that, in principle, there is enough information in the two decay rates for a clean determination 
of ImA/. The caveat is that disentangling ImA/ from 1 ± |A/p requires that all these quantities be affected by the 
same normalization. This may not be true once the experimental cuts, in particular possible cuts on the time i, are 
folded into the analysis. 

Using Eqs. we find 

± = ± i/^i' + 2E Re ± «^«;)] ' (15) 

i i<j 

where i and j can take the values c, d, and n. The notation i < j means that the (i, j) pairs are not repeated. Untagged 
decays probe the observables corresponding to the -I- sign. Tagged, time-integrated decays probe the observables with 
both signs (or, what is the same, they measure |Ap and |Ap separately). Since the functions fi are quite distinct in 
their Dalitz plot structure one can treat the coefficient of each different pair fif* as a separate observable.^ We may 
define the nine untagged observables by 



\A\' + \A\^ 



J2Uu\M' + 2j2 U^i Re (/,/;) - 2 ^ Im (/J* 



(16) 



where 



2T2 



_ Re [a^a* + a,a*) 

r Im (a,a* 4- a,a*) 
Ul,- ^ ij., (*^J)- (17) 

We will refer to these observables generically as the f7-observables. We may define similar quantities with \J replaced 
by D to describe the corresponding observables obtained for the difference \A\^ — These D-observables are 

obtained from Eqs. ( |l7|) by replacing the signs with '-' signs. 



Using Eqs. (g) and (|9|), the untagged observables are 



t/cc- ''^^^r' =ii-^p+kop, (18) 

C/..= ^^ = kP + kiP, (19) 

+ ^Ui, = ^£^I±|£^ = - ,rl - rozl (20) 

U^n + = -z*~\zf- |ro|^ (21) 



t/i + ^UL - ^^<±^ = z,rl + m*, (22) 
^»n- '""'',yJ""'' =N^ + kop. (23) 



*Note, however, that the errors in the extraction of these various observables are correlated. Thus, our observables, while 
distinct are not technically independent. A discussion of error correlations is beyond the scope of this paper, but will of course 
be important in the application of this approach to data. 
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Thus the Ui contain nine different functions of the four complex parameters Zi,ri. The requirement that the combi- 
nation Ucc + 2[/^ + Unn equals 1 yields the overall normalization T. 

Similarly, the tagged, time-integrated measurements will provide the additional observables 





Dcc^ 


\n |2 1;^ |2 

Pel — \0'c\ 




21 ^ 




Ddd^ 


1 1 9 1 — 1 9 

\ad? - \ad? _ 




21 ^ 








2r2 


Dcn^ 


- = 


flcfl* - flea* 


2r2 






a^a* — a^a*! 












2r2 



-2Rero + 2Rc(2r;;), 

2Re{zirl), 

zl - zzl - ?'or* 

rl - 2Re(zr*), 

ziz* + rirl, 

2Re{zrl). 



(24) 
(25) 
(26) 
(27) 
(28) 
(29) 



Note that Dcc+'2D^+ D„n = 0; thus only eight different combinations of T, zi, ri, z, and ro get fixed by measurements 
of the D-observables. 

Much can be learned about the interplay between the Dali tz plot analysis and CP-violation by looking at Eqs. (|T^ ) 
through (^9|). From the definitions of tq and ri in Eqs. (A3), we know that these quantities involve the product of a 
penguin contribution with sin a. Any nonzero value for tq and/or ri signals direct CP-violation. It is true that such 
quantities do not by themselves allow us to determine the size of CP- violation, because sin a appears multiplied by 
an unknown parameter; the theoretical calculation of this parameter is plagued by large hadronic uncertainties. (This 
is as expected for any observable probing direct CP- violation.) 

As expected, the quantities Dec, Ddd, -D^i j -^cn) ^^id D„„ are CP-odd. Surprisingly, the quantities [/^, J7j^, J/JJ^, 
and J7j„ are also CP-odd, despite the fact that they are obtained by looking for untagged decays. How does this 
come about? The reason is that there is a source of sensitivity to CP-violation induced in the Dalitz plot analysis by 
the fact that and p_ are CP-conjugate of each other. Said otherwise, when one performs a CP-transformation 

on /+a-| , one obtains /_a \- and not a quantity proportional to As pointed out before, this means that a CP 

transformation on ad yields —ad- Therefore, the quantities linear in ad and ad have peculiar CP-properties; C/^, C/j^, 
UL and C/i„ are CP-odd, while D^^, Di„ D^, and i^^ are CP-even. 

Additional observables are obtained in the tagged, time-dependent decays, which contain a sin Am/: term given by 



Im( ^AA* 1 



^ii I fi 



i<j 



(30) 



Here 



Re 



I^i = Im ( ^ttitt* ] /T^ 



a, a. 



(31) 



As before, i and j take the values c, d, and n, and the notation i < j means that no pair gets repeated in the 
sum. As a result, we have nine new observables. These observables depend on sin2(a -|- 9d) and cos2(a -I- 9d), with 
coefficients given in table I. 

We stress that the uncertainties associated with the exact shape chosen for the p resonances are likely to affect the 
Uii, Da, and In coefficients less than they affect the coefficients Uij, Dij, and lij with i ^ j 



III. ANALYSIS 



This section reviews what can be learned in the various experimental searches. Eight of the ten parameters needed 
to extract a can be fit with untagged decays alone. This greatly increases the data sample that will be available for 
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determining these parameters, since at the i3-factories tagging efficiencies are estimated to be of order 0.3 (or less). 
Further, many questions about backgrounds and the contributions of the other resonances to the three pion Dahtz plot 
can also begin to be answered using this larger data sample of untagged events; though they must also be re-examined 
in the tagged data sample, where non-_B background will presumably be reduced. Fitting to tagged, time-integrated 
events, fixes one further parameter and gives eight additional measurements that depend on combinations of the 
eight parameters already fixed, thus improving the precision of their determination. Only the important CP- violating 
CKM-related parameter a + 9d remains to be fit to the tagged time-dependent Dalitz plot data. Our conclusion is 
that these preliminary steps can and should be performed at both symmetric and asymmetric i3-factories. 



A. Observables from untagged decays 

The observables Ucc through [/„„ can be combined to yield T, z, ri, \zi\, |ro| and arg(2irg). Therefore, the 
nine observables present in untagged decays, Ucc through J7„n, allow us to measure eight quantities and give one 
mathematically redundant piece of information, which of course serves to further constrain that combination of 
observables. 

The result of this analysis is that the large data sample of untagged decays is extremely important for the final 
determination of a in the B ^ pn channels. One may use untagged decays to measure all relevant quantities except 
one angle — the angle between zi (or rg) and z (or ri) — and the CP- violating phase a + 9d- Moreover, one will be 
sensitive to direct CP- violation through |ro| and |ri|. 



B. Observables from tagged time-integrated decays 



The subset of events corresponding to tagged, time-integrated decays will provide measurements of the additional 
observables Dec through Since Unn and are expected to be small, it is interesting to note that the remaining 

16 U- and D-observables are sufficient to fix the nine parameters T, zi, ri, z, and rg, and give seven mathematically 
redundant pieces of information. That means that one does not need to probe quantities quadratic in |a„| and |a„| 
(that is, one does not need to have an observable p^ir'^ branching fraction) in order to determine all the observables 
attainable with these measurements. Should Unn and Dnn be measured, they will provide two further mathematically 
redundant pieces of information. All these nine parameters have model-independent information that, like a branching 
fraction, can be taken from one experiment and used as input in another. Therefore, these experiments can and should 
be performed both at CLEO and at the asymmetric P-factories. The advantage of this is that fewer parameters remain 
to be determined from the relatively small time-dependent data sample. 



C. Observables from tagged time-dependent decays 



Once a data sample of tagged, time-dependent events is available, their rates contain a sin Ami term, allowing for a 
measurement of Ice through /„„. (One expects that, for some time to come, this data sample will be small compared 
to the data samples discussed above.) The J-observables depend on sin2(Q; + 9d) and cos2(q! -|- 9d), with coefficients 
given in table I. We have already seen that the combination of untagged and tagged time-independent decays yield T, 
zi, ri, z, and tq. These must now be combined with the observables in Ice through /„„. A combination of any pair 
of these /-observables yields sin2(Q! -I- 9d) and cos2(a -I- 9^) independently. A fit to all the observables determines 
2(a + 9d) (up to discrete ambiguities) and provides additional mathematically redundant pieces of information. 

We should point out that, as expected, one cannot extract a-\-9d unless some information is known about quantities 
linear in a„ and a„. However, since these amplitudes into neutral p^n'^ might be color suppressed, it is interesting 
to ask what one may learn while only bounds, rather than measurements, on \an\ and |a„| are known. In section IV 
we will show that bounds on these quantities can be combined with measurements of quantities involving the charged 
p^n^ channels in order to determine a + 9d, up to an error that decreases with |a„| and |a„|. 



IV. TWO USEFUL BOUNDS 



In this section we suppose that the bands in the Dalitz plot corresponding to P — > p^ir^ have been measured while 
only bounds on the pieces linear and quadratic in an and a„ are known. We will show that one may still find bounds 
on a + 9d, with an error that decreases as |a„| and |a„| decrease. 
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A. One useful bound 



As we have seen before, a measurement of quantities which are independent of a„ and a„ is not enough for a 
determination oi a + 9^- In particular, we can see from Eqs. (^8|), (|2j) and the first entry in Table 1 that measuring 
the parameters Ucc, Dec, and Ice, which refer only to the decays into p^7r^, is enough to determine^ 



Im 



q Oa 

P la 



= sin [2{a + 9d + So,)], 



(32) 



where we have defined 



25^ 



arg 



1 - z + ro 



1 



ro 



(33) 



Unfortunately, da is neither zero nor calculable. 

However, as we will now show, bounds on |a„| and |a„| are enough to constrain da- As these bounds decrease, the 
deviation of Eq. (p3) from sin 2 (a + 9d) also decreases. 




FIG. 1. Triangle construction for 1 — ziro and 2±ro. These complex numbers are related to Oc, flc, a„, and a„, respectively. 



To prove this, we use Fig. 1 to derive 

\2rof = \l-z + rof + |1 - z - rof - 2|1 - z + ro| |1 - z - ro| cos (2(5„), (34) 

and 

|2ro| < |z + ro| + |z-ro|. (35) 
Using these equations and the expressions for |a„|, |a„|, jad, and jad in Eqs. ^ and (^, we find 

sin^ Sa < + (36) 
A\acac\ 

Notice that, since we assume jad and jad to be known, we only need to know bounds on quantities linear in a„ and 
a„ in order to get bounds on |a„| and |a„|. The observables quadratic in a„ and a„ are not needed. We have thus 
proved that one can determine (a + 9^) by measuring only quantities from charged p^n^ final states, up to an error 
that decreases as the bounds on the color suppressed amplitudes B — s- p^Tr" decrease. 

It is interesting to compare the bound in Eq. (^6|) , with similar bounds obtained previously for the _B — > tttt decays 
[0-^ . We note a similarity between the parametrization of the B —> tttt decays and the parametrization of the B ^ pir 
decays, related through 

a„ ^ 2A(B" ^ ttV") , 

<^ \/2A(B° ^ TT+TT") , (37) 



^Notice that , ~ I 1 + i?Sf I . Therefore, measurements of Ice and Ucc determine sin \2(a + 9d + Sa)], up to an 

error that is of second order in Dec (second order in ro). 
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and similarly for the complex conjugated channels. After these substitutions and some straightforward algebra we 
can write a bound similar to that in Eq. (^ ) as 



cos (2(57r7r) > 



(38) 



where 



adii 



B{B" TT+TT-) - B{BO TT+7T- 

B{BO TT+TT-) + B{W ^ vr+vr- 



(39) 



If there is only data on untagged B —^ vr^vr'^ decays, we may still obtain a bound by substituting the squared quantity 
on the right hand side (RHS) of Eq. (||) by 2B{B" TT°n") + 2B(B° ^ tt^tt^), a bound previously derived by Charles 
^tij,^ addition, there is no data on adir, then we may obtain a weaker bound by setting Odir to one on the RHS of 
Eq. (^) p. In this form, the bound depends only on untagged data and it is related to a bound obtained earlier 
by Grossman and Quinn Q by using B^ tt+tt" decays instead of i? tt+tt" decays on the RHS of Eq. (|38|). 



B. A bound from interference effects 



One may find a much cleaner bound by rewriting Eq. (|3^) in the form 

1 + x 



26a = arg 



1-x' 



where 



We find 



1-z U,, + UR + iUl^ 



tan {25a) 



2Ima; 



(40) 



(41) 



(42) 



If the quantities z and rg are of order one, then we expect to be able to measure them. The case of interest here is 
when z ^ 1, ro ^ 1, and the best we can do is place bounds on quantities linear in these variables, while we expect 
to be able to measure Ucc and Ice- In this case, we can constrain the allowed values for tan {25a) by combining the 
measurement of Ucc with the bounds on J7^, J7j„, and -D^„. (In effect, to gain any information, we need the 
bounds on the magnitudes of the latter quantities to be smaller than the value for Ucc-) 

Notice that Eq. (41) involves U^, J7j„, I?^ and Z3^„ which are linear in a„ and a„. That means that here one is 
using the interference between the tails of two different pir channels. This has two consequences. Firstly, this bound, 
which is extremely powerful, has no analogue for the B — )■ tttt decays. Secondly, since the bound depends on the 
interference effects, it may be sensitive to the assumptions mentioned above about the exact shape of the p-resonances. 



V. SOME EXPERIMENTAL ISSUES 



This analysis has taken a purely formal approach and has not evaluated all the relevant experimental questions. The 
quantities which we call "mathematically redundant" are actually additional data samples that contribute to fixing the 
parameters. The parameters we define are all intrinsic to the physical process and hence are, like branching fractions, 
expected to be the same in any experiment. One issue that will be important experimentally is that the errors on 
the various parameters are highly correlated and must be treated correctly in establishing bounds such as those from 
Eq. (42) and the allowed range for a. Moreover, backgrounds, efhciencies, and cuts will differ in the different data 
samples and must be investigated separately in each case. A good knowledge of the variations in efficiencies across the 
Dalitz plot is also essential for this analysis. The extraction of the parameter T is sensitive to the knowledge of overall 
efficiencies; it may be better to simply define a T^ff measured separately for each data sample than to depend on the 



9 



accuracy with which the overall efficiency and luminosity for each data sample is known. Other caveats have been 
mentioned, such as the need to explore the sensitivity of the results to reasonable changes in the assumed p-shape 
parameterization, and the recognition of the possible large backgrounds in the untagged sample. All of these issues 
and many more will only be settled by examining the data.^ None of them appear to us to invalidate the expectation 
that it will be very valuable, at least in the early years of study of these modes, to use the parameters determined 
from the time-integrated experiments when studying the CP violating effects in the time-dependent data. 



VI. CONCLUSIONS 

The B ^ pir decays are described by ten parameters, c.f. Eqs. (||) and @. We have shown that untagged data 
can be used to extract eight of these ten parameters and tagged time-integrated data allows evaluation of one further 
parameter. This leaves only the one CP- violating parameter a {a + 9d if new physics contributes an extra phase 6d to 
the mixing) to be determined from time-dependent data. The parameters in question are defined in an experiment- 
independent fashion and hence the values measured in a time-integrating experiment, such as can be pursued at 
CLEO, can be used as input to fits of the time-dependent data sample - thereby possibly allowing a measurement of 
the parameter a + 9d earlier than could be achieved without this input. 

We have also shown that, if the neutral p contributions are small, then, prior to the time when the statistics are 
sufficient to provide a measurement of these effects, bounds on their contribution will allow bounds on the shift of the 
angle measured from the rates and interference of the two charged p channels from the true value a. 
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APPENDIX: PARAMETRIZING THE DECAY AMPLITUDES 

In this appendix we parametrize the B*^ — > tt^tt^tt^ decay amplitudes by breaking them into CP-even and CP- 
odd components. This will allow us to see more clearly what quantities may be measured with the various types of 
experiments. We may decompose the isospin amplitudes into tree-level and penguin contributions as 



3/2,2 
3/2,1 



-P 



3/2,26 



27 



-P 



3/2, le 



^1/2,1 = -Pi/2,ie' 



Pi 



1/2, le 



^1/2,0 = -Pi/2,oe*'' + ^'i/2,oe 



1 



(Al) 



where we have used the fact that ^3/2,1 and ^3/2,2 only receive contributions from the tree-level diagrams. For 
convenience, we have included an explicit minus sign in the definitions of the tree-level amplitudes. Substituting into 
Eqs. (|l|) and (|^), and dropping the A/ = 5/2 amplitude, we find 







/2Pl/2,0 






3 P 


ioL 


^3/2,1 


, ^1/2,1 




P 


P 


ia 


l-v 


/2Pl/2,0 






3 P 



2Pi/2,o , . /2 ^1/2,0 . 

— cosa + zy g^^sma 



3 

^1/2,1 



P 

2Pl/2,0 

3 P 



^ .^1/2,1 
cos a + t — ^ — sm a 



T 



cos a 



2Pl/2,0 

3 P 



sm a, 



(A2) 



One can get some idea of the impact of some of them by simulations based on model values for the parameters. Such studies 
are in progress [ p^ . 
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where we have defined T = \/3T3/2,2. 
Let us define 



2 /2T1/2.0 /2P1/2.0 

^3/2,1 , ^1/2.1 , ^1/2,1 

21 = —jr- + + cosa, 
2P1/2.0 . 



■-^1/2,1 . ,f,n\ 

ri = I — ^ — sma. (Aoj 

The parameter zi (z) contain the CP-even contributions to the final state with isospin If — 1 {If — and also 
// — 2), while ri (ro) contains the CP-odd penguin contributions to the final state with isospin If — 1 {If — 0). 
Using these definitions, we arrive at 

e^l'a, = Te-*" (1 - z - ro) , 
e'% = re-*" (zi+n), 

e'f'a^^Te-'^iz + ro). (A4) 

Eqs. (||) are obtained from these by dropping the (irrelevant) overall phase factor e"^^ . Similarly, we may reach Eqs. (||) 
by applying CP-conjugation, multiplying by q/p, and removing the same overall phase factor e'^^ . 

We may also parametrize the amplitudes involved in the — > p°7r^ 7r+7r~7r^ and p*7r° tt^tt^tt'^ 

decay chains as 



e 



'^AB±-,pO^± = — j=- I - - zi - SiTri] , (A5) 



and 



V2 V2 



e^"^AB±-.p±^o ^ i^- + zi+Si±nj , (A6) 

respectively, where we have dropped the A/ = 5/2 amplitudes. These decays involve the new complex parameter 

^1 = -^^. (A7) 

Therefore, the information in B^ — > p^tt^ tt+tt^tt^ decays by itself does not help in constraining the parameters 
involved in the decays of the neutral P's into three pions; these two rates merely allow a determination of the 
magnitude and phase of 5i . The decays from charged B mesons are only useful if one can measure both the B^ — > 
p°7r^ 7r+7r~7r^ and the experimentally challenging B^ p^n^ Tr^Tr^Tr" decay chains. For example, one might 
use the information from the decays of the neutral B to get T, ri, and zi, and the B^ — > p'^n^ decays to get ^i. One 
would then be able to predict the rates for the B^ p^n^ decays. Even with the additional direct CP asymmetry 
measurements in these channels, we see that we have no way of extracting a measurement of sin a. As before, only 
the quantity ri appears. 
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Observable 


Coefficient of 


Coefficient of 




sin 2(a -he'd) 


cos2(q + Od) 


Ice 




2Im [ro(l-z)*] 


Idd 




2Im [rizl] 


led 


2Re[ri{l - zy +zir^,] 


2lm[{l~z)zt+rort] 


-'cd 


2Im [ri(l-z)* +2iro] 


2Re [(1 - z)zt + ror^] 


/J 

cn 


2Re [z - \z\'^ + Irol^J 


2Im[(l - 2z)ro] 


J« 

cn 


21m z 


2Rero 


ji 

'-dn 


2Re [zri — rozl] 


2Im [rirg — z\z*\ 


'dn 


2Im [zrl ~ rozl] 


2Re [rirS — z\z*\ 


Inn 




2Im \ztq\ 



TABLE L 7-observables obtained from the sin Am t term. The table shows their dependence on sin2(a + ^d) and 
cos2(a + ed). 
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